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The possibility that like-charges can attract each other under the mediation of mobile counterions 
is by now well documented experimentally, numerically, and analytically. Yet, obtaining exact 
results is in general impossible, or restricted to some limiting cases. We work out here in detail 
a one dimensional model that retains the essence of the phenomena present in higher dimensional 
systems. The partition function is obtained explicitly, from which a wealth of relevant quantities 
follow, such as the effective force between the charges or the counterion profile in their vicinity. 

Isobaric and canonical ensembles are distinguished. The case of two equal charges screened by an 
arbitrary number N of counterions is first studied, before the more general asymmetric situation is 
addressed. It is shown that the parity of N plays a key role in the long range physics. 


I. INTRODUCTION 

Coulombic effects are often paramount in soft matter systems, where the large dielectric constant of the solvent 
(say water) invites ionizable groups at the surface of macromolecules to dissociate [THS]. While a realistic treatment 
requires considering three dimensional systems, interesting progress has been achieved for lower dimensional problems 
where the key mechanisms can be studied in greater analytical detail [H]. In particular, a one dimensional model 
was introduced in the 1960s by Lenard and Prager independently, for which a complete thermodynamic solution was 
provided [319]. This model has been further studied in Ref. m, but it turns out that some interesting features have 
been overlooked in relation with the like-charge attraction phenomenon min]- This striking non mean-field effect, 
relevant for strongly coupled charged matter mma is the thread in our study. 

The paper is organized as follows. The model is first defined in section [Ilj It mimics the screening of charged 
colloids. The Coulomb potential in one dimension between two charges q and q' located along a line with coordinates 
X and x' is 


v{x,x') =—qq'\x — x'\. (1) 

Therefore, the electric field created by one particle is of constant magnitude. This fact simplifies the study of the 
equilibrium statistical mechanics of such systems, and allows to obtain some of its properties by simple arguments. 
Furthermore, it also allows for an explicit computation of the partition function miH]. The system under scrutiny can 
be envisioned as a collection of parallel charged plates, able to move along a perpendicular axis. The salient properties 
of this system can be obtained by simple arguments which we present in section [Ilj followed afterwards by a more 
technical analysis where the explicit calculation of the partition function is performed, first in the isobaric and then 
in the canonical ensemble. After having presented the symmetric case, section m will generalize the investigation 
to the situations where the two screened charges are different. Noteworthy is that parity of the particle number 
considerations will play an important role in the remainder. 


II. SCREENING OF TWO EQUAL CHARGES BY COUNTERIONS ONLY 

Consider two charges q along a line located at Y = 0 and x = L. Between the charges there are N counterions 
of charge e = —2qfN between them. Consider the equilibrium thermal properties of this system at a temperature 
T, and as usual define P = l/{kBT) with ks the Boltzmann constant. This simple model mimics the screening and 
effective interaction between two charged colloids in a counterion solution, without added salt. In one dimension, Pe^ 
has dimensions of inverse length, therefore it is convenient to use rescaled units in which all distances are measured 
in units of l/(/3e^): x = Pe^x. It is also convenient to work with a dimensionless pressure P = P/e^ where P is the 
pressure (equal to the force, in one dimensional systems). 

The potential energy (dimensionless, measured in units of ksT) of the system is 
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FIG. 1: An odd number of mobile counterions screening two like-charges. The N mobile ions (counter-ions) have charge —2q/N 
and the conhning objects have charge q, so that the whole system is electro-neutral. Here, N = 2p-\-l is odd, so that a single ion 
(referred to as the misht since the net electric force acting on it vanishes) “boats” in between the two screened boundaries which 
attract, each, p ions in their vicinity (see also Fig. |^. This single free counterion provides the binding mechanism responsible 
for long range attraction. In the canonical treatment, L is held hxed, while in the isobaric situation, it is a fluctuating quantity. 
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FIG. 2: An even number of counterions screening two like-charges {N = 2p). At large distance, the two double-layers (made up 
of an ion q and p counter-ions) decouple since they are neutral. No misht ion is present to mediate attraction, and the pressure 
is repulsive at ah distances. 


Before presenting the technical analysis, we start by simple and more quantitative considerations. 


A. Possibility of attraction between like-charges 

1. A heuristic argument 

The possibility of attraction between the two -^q charges at 0 and L is related to the parity of A. If A is odd, 
A = 2p + 1, then p counterions will form a double layer around each charge q. This will form two compound objects 
with charge q{l — 2p/N) = q/N each one, located around 0 and L. There will be in addition one counterion between 
these two object, which is essentially free, as the electric field created by the charges located on each side around 0 
and L cancel each other. When L is large enough, consider figure The right side of the system composed of one 
charge q and p counterions has charge q/N. The left side which, for the sake of the argument, has the free counterion 
plus the compound charge, exhibits a total charge —q/N. Thus the force exerted by the left side on the right side is 
P -q^/N^ = -eV4, an attractive force. Thus one expects that P ^ —1/4, for L ^ oo. 

On the other hand, if N is even, there will not be a free counterion between the layers, which will be completely 
neutral, thus one expects that P ^ 0+ when P ^ oo, as shown in figure 


2. Beyond heuristics 

The previous intuition, providing a large distance attraction for odd A, can be substantiated by a simple calculation. 
Use will be made here of the contact theorem [iniiisHis], an exact relation between the force exerted on the charge 
g, and the ionic density at contact (stemming from the mobile charges —2q/N). Such a relation is particularly useful 
for discussing the like-charge attraction phenomenon nainiiii]. The argument allowing to get the contact density 
is two-fold, and goes as follows. 
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FIG. 3: Upon regrouping the p+ 1 leftmost counterions in Fig. one obtains an ion with charge —q — q/N. This newly defined 
system has the same large distance pressure as that of Fig. 


First, we argue that at large L, the p counterions that are closest to each boundary remain in their vicinity, while 
the middle free counterion (the misfit in Figs. and , which does not feel any electric field by symmetry, tends to be 
unbounded and no longer contributes to the pressure (discarding 1/L terms). In a second step, we thus compute the 
contact density in a system of an isolated charge with a double-layer of p ions in the vicinity (the total charge of 
this composite object, shown on the right hand-side of Fig. is q/N. The solution to this problem is not immediate, 
but can be found by a convenient mapping onto a more convenient problem, shown in Fig. As illustrated in the 
figure, we regroup the p + 1 leftmost counterions in a single ion, having charge —g(l + '^/N). At large distances, 
this regroupment does not influence the distribution of counterions around the rightmost ion and thus leaves the 
large L pressure unaffected. The next important argument is that the pressure can be equivalently computed from 
the contact density at the rightmost, or leftmost charge -\-q. It is thus simpler to perform the calculation in the newly 
defined regrouped system (left hand side of Fig. |^. The regrouped ion with charge —g(l + '^/N) is in the electric 
field of the charge q on its left, and of the composite system on its right having charge q/N. This amounts to a field 
^'(l — l/N). Hence, the electric potential energy reads g^(l — 1/N)x{l + l/A"). The corresponding Boltzmann weight 
gives the density of the regrouped ion 


p{x) 



exp 




( 3 ) 


where due account was taken of normalization (f pdx = 1). The contact density p(0) = f3q‘^{l — 1/N^) finally yields 
the pressure through the contact theorem /3P = p(0) — (dq^. We get here P = —(^/N^ (or equivalently P = —1/4), a 
result which by construction holds in the large L limit. The reason for a non vanishing pressure at large distance is 
that the p counter-ions cannot exactly screen the charge of an ion q. It is no longer the case when N is even, in which 
case P ^ 0 for L ^ oo. The present results will be fully corroborated by direct partition function calculations. 


3. Correction to large distance asymptotics and crossover pressure 


Returning to the case when N = 2p-\-l is odd, we can also estimate the first correction to the pressure for large L. 
Consider that L is fixed (canonical ensemble) and large. Since the system is somehow equivalent to two double layers 
with a free counterion in between, this counterion will contribute to the pressure (denoted as Pc in the canonical, 
fixed-!/ ensemble) with a correction 1/L. This estimate can be made more quantitative. The available space for the 
free counterion is not L, but it is rather L minus the space occupied by the diffuse counterion layers, given by {xp)oo the 
thermal average position of the p-th counterion if they have been ordered Xi < X 2 < ''' < Xp < < • • • < X 2 p+i, 

in the limit L ^ oc. Thus 


4 + L-2(rCp)oo ■ 

This is illustrated in figure]^ In the following section, we evaluate explicitly {xp)ao and find 
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FIG. 4: An odd number of counterions screening two like-charges. The free “misht” ion is singled out. 



FIG. 5: The (canonical) pressure Pc as a function of the separation L, for A" = 25 (continuous bottom line) and A = 26 
(dashed upper line). For A odd the pressure becomes negative at large distances. 


In the other limiting case L ^ 0, the result is [TO] Pc = N/L, that can be understood as all the A counterions are 
squeezed in a small distance L. Thus we see that the pressure is positive (repulsive force) for small separations L ^ 0 
then changes to negative pressure (attractive force) for large L. 

We will show in the following section that the o{l/L) corrections in (§ are actually exponentially small, in the 
canonical ensemble, therefore equation §) gives a fairly good approximation for the pressure for a large set of values 
of the separation L. From this, one can estimate the distance L*, at which the effective force between the two charges 
becomes attractive 


A - 1 

L* c^A + 2{xp)o.=4 + 2j^. (7) 

Figure shows the pressure Pc as a function of L, for A = 25 and for A = 26 particles. For A = 25 (odd) the 
pressure changes its sign at L* = 4 + 2 * 24/26 ^ 5.85, while for A = 26 the pressure is always positive. 

Summarizing, in the case of odd A, the possibility of having an effective attraction for large separations L is due 
to the sharing of the “free” ion which leads to the creation of opposite charges objects (ions q plus their counterion 
clouds). Although the analytical results presented here are valid only for this one-dimensional model, the same 
physical mechanism has also been observed in three dimensional systems [191 [2D] • It can also be surmised that in 
situation of odd A where the free counterion has a varying charge, attraction will be all the stronger as the charge 
will increase in absolute value. In addition, the very mechanism brought to the fore here indicates that at mean-field 
level, where the discrete nature of ions is discarded, attraction should be suppressed, which indeed is the case EIHIS]- 
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B. Explicit exact calculation of the partition function 

1. Preliminary observations 


The equilibrium thermodynamics of the one-dimensional two-component Coulomb gas was solved simultaneously 
but independently by Lenard [7] and Prager [8]. In the present model, only one type of identical particles (the 
counterions) are present. It is convenient to order the particles as 0 < < • • • < < T. Then, rearranging the 

terms in the potential energy of the system can be written as 


and 



p-i 

2 '^ip - j){x 2 p+i-j - Xi+j) for N 
j=o 


2p + 1 odd. 


N‘^L ^ ^ 

U = — -y](2p - 2j - l){x 2 p-j - Xi+j) 

j=o 


for N = 2p even. 


(8) 


(9) 


Notice that in the case N = 2p + 1, the particle with position does not appear in the potential energy. It is 
the free counterion (misfit) discussed in the previous section, whose role is crucial for the possibility of like-charge 
attraction. 

The canonical configuration integral is 


pL pXN rX3 pX2 

Zc{N,L)= / dxN / • • • / dx2 / dxie~^ . 

Jo Jo Jo Jo 


( 10 ) 


As mentioned by Lenard in his seminal paper [7] “the (configuration) integral is elementary (because) the class of 
functions consisting of exponential of linear functions is closed under the operation of indefinite integral (...) however 
the task of evaluating (it) is not trivial”. For small N one can compute by hand Zc, and for larger given values of 
N it can be obtained numerically with the aid of a computer algebra system software program. By inspection of the 
integral (10), one can deduce that Zc is a linear combination of products of exponentials of L and linear functions of 
L. One can also deduce the argument of each exponential function of L by keeping track of the factor that multiplies 


each Xk in the integral (10). These come from the explicit term in U (for instance, for it is 2{p — j) in the case N 
odd), but after each successive integration, the factor of Xk will by added to the one of Xk+i due to the upper limit of 
integration. Taking that into account, one realizes that the exponentials of L in Zc are of the form exp( —(j + |)^A) 
in the case N odd, and exp(—j^L) in the case N even. Thus, the configuration canonical integral is expected to be 
of the form 


and 


p 

Zc{N,L) = ^e-^^+^^^^{AjL +Bj) for AT = 2p + 1 odd, (11) 

i=o 


P 

Ze{N, L) = J2 + Dj) for N = 2p even. (12) 

i=o 


The non trivial task is to evaluate explicitly the coefficients A 


3^ 


Bj, Cj and Dj. 


This is done in section IIB 5 


2. Previous results 

In [To], the present system was studied, but an exact analytical explicit evaluation of the partition function for an 
arbitrary number of particles was not achieved. Rather, an interesting reformulation of this model was proposed, by 
mapping it into a quantum mechanical problem, following a technique put forward by Edwards and Lenard [9|. It 
was shown in m that the configuration integral is given by 


Zc{N,L)=b{N/2,N/2,L) 


( 13 ) 
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where b{n, N/2,x) is the solution of a set of N coupled elementary linear differential equations 

db{n, N/2, x) 


dx 


= —(n^/2) b{n, N/2^ x) + b{n — 1, N/2^ x) 


with the initial condition 6(n, A^/2,0) = Sn-Nl 2 - Integrating this equation one has 

b{n,N/2,Xn)= [ - l,N/2,Xn-i) dxn-i. 

Jo 


(14) 


(15) 


Then, starting from the known b{—N/2^ N/2, xi) one has to perform successively N integrals (15) to obtain 
b{N/2, N/2, L) and the configuration integral. This task is equivalent to performing directly the N integrals of 
the configuration integral (10). Thus, unfortunately, the method proposed in m does not provide any computational 
advantage over a direct numerical evaluation of the partition function. 

Here, our goal is to obtain an explicit analytical expression for the configuration integral for an arbitrary number 
of particles N. Using Lenard [7] and Prager [8] method, we will first compute the partition function of the constant 
pressure ensemble 


Zp{N,P)= e-^^ZciN,L)dL (16) 

Jo 

which is the Laplace transform of the canonical configuration integral Zc. This is a straightforward application of 
the technique of Lenard and Prager, and it is actually much simpler than the complete work presented in [71 [8], 
since all particles are identical and we will not have to deal with the combinatorial problem of studying the different 
configurations of charges. 

Then, we shall invert the Laplace transform to obtain the canonical, constant “volume” L, configuration integral 
Zc{N^ L). Since we are interested in finite systems, the results from the canonical ensemble and the constant pressure 
ensemble will differ, and it is of interest to compare them. 


3. Evaluation of the dijfuse layer size {xp)oo 

To introduce the technique used to compute the partition function, we undertake in this section a preliminary, 
simpler task, based on the same technique: the exact evaluation of the diffuse layer size {xp)oo- This quantity 
appeared in the discussion of section [IlA[ Consider here that L ^ oc and N = 2p+ 1. The double layer composed by 
the charge q at L and its corresponding p counterions are thereby ‘sent to infinity’. The remaining p + 1 counterions, 
however, still feel the electric field created by this far charged double layer. The potential energy part which depends 
on the position of the remaining counterions is 

p-i 

Ucc = - j)xi+j . (17) 

j=0 


We wish to evaluate 


^ U.l=ldXk 
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(18) 


Let 

F{s)= f e-^^-^^^/‘^dxi...dXp. (19) 

J 0<X\<X2<---<Xp 

Then {xp)oo = —2dlnF{s)/ds\s=o- Following Lenard [7] and Prager |8] it is convenient to re-write the potential 
energy as 


U^ = 


X] ((p - J + 1)^ + {p-j + 2)^) {xj - Xj_i) - : 


1=1 


( 20 ) 
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with the convention that a;o = 0. Let us define 

fj{x) = e-[(P-J+i)"+(P-J+2)"]^/2ii-(a;) 
where H{x) is the Heaviside step function. Then 

roc P 


POO POO Jr' 

F{s) = / dxi--- dXpY\fj{xj-Xj-i)e 
Jo Jo 


— {s — l)xp/2 


( 21 ) 


( 22 ) 


We notice that F{s) is the Laplace transform (evaluated at {s —1)/2) of the p-fold convolution product /i * /2 * * * * * /p- 
The Laplace transform Cfj of fj is elementary 
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(p - j + 1)2 + (p - j + 2)2 + S - 1 2{p- j + l)(p - i + 2) + , 


Then 
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Computing the derivative of lnF(5) we obtain 

dlnF(s)| ^ 
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ds 


^ 1 
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*=o ^ (p-i + i)(p-i + 2) 


y _ I _= 
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Thus proving 


pFlJ p +1 


4- Isobaric ensemble 


(23) 


(24) 


(25) 


Consider now the finite system with L < oc. We will detail the calculations in the case = 2p + 1 odd, the case 
N even can be obtained by a simple adaptation of the same technique. As it was done in the previous section, it is 
convenient to re-write the potential energy ® as 


L 1 


p+i 


u = -^ + {{p - j + y + {p - j + 2 )^) {x2p-j+3 - X2p-j+2 + Xj - Xj-i) 


(26) 


i=i 


where, by convention, we defined xq = 0 and X 2 p -\-2 = L. With fj defined in (21), we notice again that the canonical 
partition function is a convolution product of 2p + 2 functions fj 


Z,{2p + 1, L) = fj * ffj (L) . 


(27) 


The isobaric partition function Zp is the Laplace transform of Zc, and we have 

p+y . . N X 2 

{ cj A p - ^ 


p+i / / . 

zp{2p+i,p) = WUfAp-j 
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where s = (4P + l)/2. Factoring (k + + P = (k + ^ — i^/P){k + 5 + i^P) = \k + ^ + i^/P\^, the above product 

can be expressed in terms of Gamma functions 


( ' V 

r(| + i/p) 


T{p+l+iy/P) 


Zp(2p + l,P) 

The average length of the system is given by the usual thermodynamic relation 

dlnZp 2 ^ 2 


(L) = - 


dP 
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(29) 

(30) 

(31) 


where 2 p{z) = dlnr{z)/dz. We can notice that this expression has a pole for P = —1/4, from which we obtain the 
behavior when (L) oc, P —1/4, in agreement with the general discussion of section [II A When N is even this 
pole is absent (see below). 

If N = 2p is even, similar calculations lead to 


Z,{2p,L) = e^/Vp+l * (.!*/i+l * fj+h) 


and 


Zp{2p,P) = ^f[ 


p 1^1 (fc 2 + p )2 


1 

P 


T{l+iVP) 


r{p + 1 + iy/P) 


(32) 


(33) 


Notice an important difference in the analytic structure of the partition function in the case N odd (27)-(28) and N 


even (32)-(33): for N even, there is a single function /p+ 3/2 in the convolution product, leading to a pole of order one 
for P = 0, m contrast to the case N odd, where the functions /p+i appear twice in the convolution product and the 
pole for the smallest value of \P\ is of order two and it is for P = —1/4, rather than P = 0. In the case N even, the 
term /p+i * /p+i corresponds to the coupling of the left diffuse layer with the free counterion and the coupling of this 
same free counterion with the right diffuse layer. On the other hand in the case N odd, the term fp-^sj 2 corresponds 
to the direct coupling of the left and right diffuse layers. 

The average length, for N = 2p even, is 


{L) 



k=l 


2 

pTp ' 


(34) 


We note that (P) ^ 00 when P ^ 0+, in contrast to what happens when N is odd, where (P) ^ 00 when P ^ —1/4. 


5. Canonical ensemble 


We return to the case = 2 p + 1 odd. To compute the canonical partition function, we need to invert the Laplace 
transform computed in the previous section 




Z^{2p+l,L)=C 
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n 
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This rather technical part of the analysis is presented in Appendix]^ where it is shown that 
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From this expression, we obtain the canonical pressure Pc = , 
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For N = 2p even, the results are 
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(37) 


(38) 


(39) 


6. Limiting cases and comparison between the ensembles 


With the exact expressions obtained above, we can prove rigorously the limiting behavior of the pressure when 
L ^ oc and L ^ 0 discussed in section nu 

Let us consider first the case = 2 p + 1 odd. In the canonical ensemble, the beha vior of the pressure Pc when 
L ^ oc, is obtained from the term j = 0 of (36), confirming the prediction <11 of section [IlA| Furthermore, we realize 
that the next to next to leading order correction is exponentially small 
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-2 
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In contrast, when N = 2p, the pressure tends to 0 exponentially fast when L ^ oo 


Pr = 


4p^e 


2 ^-L 


(p+ + 
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2p+l 

P{P + 1) 




(41) 


The behavior of the pressure is different in the isobaric ensemble. Consider again first the case N = 2p 1. From 
, we already know that when P = —1/4, (L) ^ oo. Denoting 5 = (4P + l)/2, one can expand ([30l) for small s 


(42) 


and invert the relation to obtain P as a function of (L) when (L) -+ oo. For instance, to order 0(s), Eq. (30) is 

S p i. 


where 




(43) 


with 'HP = Y%=i ^ ^ the harmonic numbers. Inverting that relation, up to order 0({L) ^), gives 


"^=-4 + 


8S(p) 


(L) - 2 ^ ((L) - 2 ^)3 




(44) 
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Notice a factor 2 of difference in the next to leading order correction (the 0{{L)~^) term) in the pressure in the 
isobaric ensemble and the canonical ensemble. Furthermore, in the isobaric ensemble the next to next to leading 
order corrections are algebraic and not exponential as in the canonical ensemble. 

For N = 2p, the behavior of the pressure, in the isobaric ensemble, when (L) ^ oc, is 


P = 


2n 


(4) 


{L) 




m 




+ 0{{L) 




(45) 


Notice again the different behavior with respect to the canonical ensemble. Here in the isobaric ensemble, the pressure 
vanishes as 1/(F), whereas in the canonical ensemble it vanishes exponentially fast, as e~^. 

Let us study the other limiting behavior of the pressure, for small separations L. Let us focus on the case N = 2p + l 
first. It is not completely straightforward to obtain the behavior of the pressure in the canonical ensemble when L ^ 0 
directly from expression (37). Rather, it is better to return to (27), and notice that if L ^ 0, then the convolution 
product fj * fj behaves as 


fj * fj{x) = xH{x) + O(x^) 


(46) 


which is independent of j. Then, 


fj * fj 

i=i 


(x) 


„2p+l 


and 


(2p + l)! 


N 


+ 0(a;2p+2) 


Z,{2p + 1,L) = ^+0{L^+^). 


We deduce that the pressure behaves as 


L 


when L ^ 0, 


(47) 


(48) 


(49) 


a result already noticed in m- Eq. (E^ also holds when N = 2p. 

In the isobaric ensemble, when = 2p + 1, if (L) ^ 0, then, necessarily, s = (4P + l)/2 
that equation to order 0(5“^), one obtains 


oo in (30). Expanding 


P = 


+ l N{N^2) 
~{L) 12 


+ 0((P)) when (P) ^ 0 . 


(50) 


This result also holds true for N = 2p. Notice again the difference between the canonical (49) and isobaric ensemble 
( [5Q| ), where the leading term changes from N/L to {N + 1)/P. 

When N = 2p 1 is odd, the pressure changes of sign when L varies. It is positive for P 0 and negative for 
P ^ oo. We already obtained an approximation of the value P* of P when this occurs in the canonical ensemble, see 
up to exponentially small corrections. In the isobaric ensemble, one just has to put P = 0 in ( [M| ) to obtain the 
exact value 


( i *)=8 l + £ 


^ ( 2 fc + 1)2 


= 7r2 -2V;'(p + 3/2). 


(51) 


Eor this quantity, the predictions from the canonical ensemble § and the isobaric ensemble ( [^ are again different. 

Eigurej^ shows the pressure as a function of the separation, for N = 15, in the isobaric ensemble and the canonical 
ensemble. Notice that the pressure from the canonical ensemble is smaller that the one in the isobaric ensemble for 
the same separation. Eigurej^ shows the value of P* for which the pressure changes of sign as a function of N, when 
N is odd, in both ensembles. Notice again that in the canonical ensemble, the change of sign of the pressure occurs 
for smaller values P* of the separation than in the isobaric ensemble. 


III. SCREENING OF TWO UNEQUAL CHARGES 

In this section we consider a generalization of the previous model, where the two charges located at x = 0 and 
at X = P are qi and q 2 , respectively, which can be eventually different. The overall system should be neutral. 
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FIG. 6: The pressure P as a function of the separation L, for N = 15. The top continuous line represents the result from the 
isobaric ensemble, and the dotted bottom line those from the canonical ensemble. 


N 


FIG. 7: The value of the separation L* for which the pressure vanishes and changes sign as a function of N for N odd. The 
filled squares represent the results from the isobar ic ensemble, and the filled disks, their canonical counterpart. 


therefore gi + g 2 = —Ne^ e being charge of one counterion. It is convenient to introduce the notation Qi and Q 2 
such that qi = —eQi and q 2 = —eQ 2 . The electroneutrality relation is Qi Q 2 = N. The charge asymmetry can be 
characterized by the quantity a = Qi — Q 25 which allows to write Qi = {N + a)/2 and Q 2 = {N — a)/2. The potential 
energy of the system is now 

N 

U{N,L,Qi,Q2) = - \xi - Xj\ + a'^Xi + {Q2)^L . (52) 

l<i<i<iV i=l 

The overall effect of the charge asymmetry is to introduce a global electric field proportional to a (the term in Xi). 


A. Isobaric ensemble 


Adapting the ideas of section IIB 4 to the present case, we can obtain the isobaric partition function. Once again, 
the results differ depending on the parity of the number of counterions N. For A = 2p + 1 odd. 


1 


while for N = 2p even. 




1 ^ 

Zp{2p,P,Qi,Q2) = ——^- II 


(f)' + P tA [(k - !)' + P] [{k + 1)2 + P] ■ 


(53) 


(54) 




























12 


The above formulas highlight the difference between the two cases, depending on the parity of N. However both 
formulas can be summarized in a single one as 


N 


Zp(7V,P,Qi,Q2) = n 


1 


N 

n 


1 


+ p 11 {i - Q<)2 + p 

1 


n 




y^ + P 


where we defined 


Q< = ^ = min(Qi,Q 2 ) and ^ = max(Qi,Q 2 ) ■ 


(55) 


(56) 


Taking the derivative of (55) with respect to P, we obtain the relation between the average length (L) of the system 
and the pressure P in the isobaric ensemble 


N 




(57) 


If a ^ Z is not an integer (gi and q 2 are not integer multiples of —e/2), or \a\ > N {qi and q 2 have opposite signs), 
then Zp has simple poles. But when a G Z is an integer and \a\ < N, the partition function Zp turns out to have 
some double poles. This corresponds to the case when 2 Qi and 2Q2 are both positive integers. In that case it is best 
to reorder the products in ( [5^ to make those double poles more apparent. The result depends on the parity of 2 Qi 
and 2Q2 (both have the same parity). If 2 Qi and 2Q2 are odd, then Qi and Q2 are half integers: Qi = [QiJ + | and 
Q2 = LQ 2 J + The notation [xj denotes the floor function of x (largest integer less or equal than x). The isobaric 
partition function (55) becomes 


LQ<J 

ZpiN,P,Qj,Q2)= n 


LQ>J 

n 


e=o [i^+h)^ + P] «=LQ<j+d^'’" 2)^ 


and the corresponding equation of state is 

LQ<J 

w = E 




LQ>J 

E 

LQ<J +1 


1^2 . 


1 


£=0 V- ' 2> 

When Qi and Q 2 are positive integers, these expressions become 

Q< 

p n /^2, p)2 Li p+p’ 

and 

Q< o 


^ Q< 

Zp{N,P,Qi,Q2) = -'[[ 


Q> 

n 


w = 4 +e 


«=i 


£2 + P 


E 

i—Q < +1 


1 


fi + p 


(58) 


(59) 


(60) 


(61) 


B. Canonical ensemble: the partition function 


To compute the canonical partition function, one has to perform the inverse Laplace transform of the expressions 
obtained in the last section. From the above discussion, it is clear that the results will have a different analytical 
structure depending on whether the isobaric partition function has simple or double poles, that is, depending on 
whether a is an integer or not. If a is not an integer, or \a\ > N, all poles of Zp are simple poles, and we obtain 
from (55): 


-XN.L.Q,M = y-l)e < i!(W-,),r(i + H + l) 

N 

3=0 


2(j-Q<)r(j-2Q<) 


j!(7V-i)!r(j + |Qi-Q 2 | + l) 


( 62 ) 
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This formula is valid whenever 2Qi and 2 Q 2 are not integers, or if Qi and Q 2 have opposite signs ((3< < 0 and 

Q> >0). 

If a is an integer, with \a\ < N, then using (58), we obtain, when Qi and Q 2 are half integers. 


LQ<J 

Z,iN,L,Q^,Q2)= Y, 


j=0 


(LQiJ + 1 + j)KlQi\ ~ +1 + j)K[Q2\ — JY 


L- 


1 


2j + l 


i^{[Qi\ - j) -i^{\Qi\ +1 + j) + '0(LQ2j - j) -'0(LQ2j +1 + j) + 2jTiy 


LQ>J 

■ E 

3— LQ< J +1 


e-(i+ 5 )"i/ 4 (j _ - 1 )! ( 2 j + 1 )(- 1 ) 4 -LQ<J-i 

(Le<J +1 + j)KLQ>J +1 + i)KLQ>J “i)- 


(63) 


and when Qi and Q 2 are integers. 




Z,iN,L,Q„Q2) Y (Q,+j)!(Q, -j-)!(Q2+j)!(Q2 - j)\ 

L — + 1 — j) — ^{Qi + 1 + j) + '0(Q2 + 1 ~ j) ~ '0(^2 + 1 + i) + V 


2 j 


J/i 




j—Q <+1 


(Q< + j)KQ> j)KQ> i)* (Qi!Q20 


2 • 


(64) 


The two previous results (63) and (64) show the different analytical structure of the two cases which depend on the 
parity of 2Qi and 2Q2^ in particular the existence of a term independent of L in the case 2Qi and 2 Q 2 even, and 
the form of the argument of the exponentials ^ (for 2Qi even), as opposed to ^ (for 2Qi odd). However, 

both results (63) and (64) can be subsumed in a single formula as follows. Let us define 


Aj{N, L, Qi, Q2) = 


(2(g<-j))2(-i)2Q>+i 

{2Q^-j)\j\{N-jy.{\Q,-Q2\+j)\ 

'0(i+i) ' 0 ( 2 Q<—i+i) T '^(i+IQi—(^ 21 + 1 ) '0(-^—i+i) T q^-j 


(65) 


L- 


m<-j) 


for jf 7 ^ Q<, and, when ( 5 < is an integer, define 

Aq^{N,L,Qi,Q2) = 

Also, let 

Dj{N, Qi, Q2) = 


\2 • 


( 0 i!Q 2 !)^ 
j!2(j + g< + l)(-l)^+2'3>+i 


(2Q< + j + 1)!(A^ + j + l)Kiei “ Q 2 I — j — 1)! 


Then, both results (63) and (64) are equivalent to 

LQ<J IQ 1 -Q 2 I -1 

j=0 j=0 


( 66 ) 


(67) 


( 68 ) 


C. Canonical ensemble: asymptotic behavior of the pressure 


For small separations L, the results (49), Pc ^ N/L (canonical) and (50), P ^ {N P 1)/(L) (isobaric), still 
hold independently of the charge asymmetry a. Thus, the effective interaction is always repulsive at short distance, 
irrespective of the charges qi and q 2 , even in the case where these charges are opposite. Indeed, the pressure is 
dominated here by the entropy cost for confining the ions in a narrow domain. 
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The behavior for large separations L will depend on whether the charges qi and q 2 are multiples of e or not, and 
their relative signs. There are four cases to consider. 

Opposite charges. First, suppose that qiq 2 < 0, the charges at the edges have opposite signs. This corresponds 
to the case \a\ > N, and the canonical partition function is obtained with Eq. (62). From that expression, we deduce 
that for L large, the leading order is given by the term j = 0 of that sum. Therefore, the effective force is attractive 
and given by 


-(Q<)^ 


OO , 


(69) 


where here (5< = {N — |a|)/2 < 0 corresponds to the charge of the edge particle which has the same sign as the 
small ions. This result can actually be obtained by simple arguments. The small ions will be repelled by the particle 
with charge corresponding to (5< and attracted to the other edge where there is a particle with charge —e(3> with 
(5> = {N + |tt|)/2 > 0. By electroneutrality, the charge of the compound object formed by the small ions and —e(5> 
will be e(5<. The effective force between this object and the other opposite charge —e(5< is repulsive, equal to 
— (e(5<)^, thus recovering (69). Application of the contact theorem of course yields the same result, since the density 
of counterions vanishes at contact with (5< (a similar effect was reported in [24l 

Like-charges that are not integer multiples of — e. To discuss this situation, we keep in mind that Qi > 0 
and Q 2 > 0 are not integers. The small ions of charge e will be divided into two parts that will try to screen the 
charges qi and g 2 - A number [QiJ of counterions will partially screen qi and [(^ 2 ] ions will partially screen the 
other charge g 2 - Each edge, with its screening cloud of counterions, will have a charge —e{Qi — [QiJ) = —e{Qi} and 
—e{Q 2 — LQ 2 J) = —^{Q 2 } respectively, where {x} := x— [x\ denotes the fractional part of x. However, since Qi and Q 2 
are not integers, we have [QiJ + [Q 2 J = N — 1: there is still one counterion to take into consideration. This counterion 
feels the electric field created by the charge difference —e{{Qi} — {Q 2 }), therefore it will be attracted to the edge 
which has the largest remaining charge (in the sense of the largest between {Qi} and {Q 2 })- To fix the ideas suppose 
{Qi} > {Q 2 }‘ The remaining ion will become part of the screening cloud of ^ 1 , and the charge of that compound 
object will be —e({(5i} — 1). Then the effective force between the two edges will be e^dQi} — 1){Q2} = —e^{Q 2}^5 
the last equality coming from the fact that {Qi} + {Q 2 } = 1- Summarizing, in general we expect an attractive force 
at large separations given by 


-(inin({(3i},{Q2})) 


OC . 


(70) 


This can be verified by identifying the largest argument of the exponentials in the canonical partition function 
or, equivalently, the largest pole of the isobaric partition function (55). The poles of the isobaric partition function 
are —{i — Q<)^, with £ going from 0 to N. Then, one can notice that £ — (5< varies from —(3< < 0 up to (5> > 0 
by integer steps of 1. From this one-dimensional array of points, we are interested in the one that is the closest to 
0. That is precisely min ({Qi}, {( 52 }), in agreement with ( [^ . One can also notice from (62) that in the canonical 
ensemble, the next to leading order correction to (70) is exponentially small of order 0 (e“IWi}-{Q 2 }|i/^ 

Like-charges that are half-integer multiples of —e. A degenerate case of the previous situation is when Qi 

In this case the canonical partition function is given by (63) 
-1/4. But the correction to leading order is 


and Q2 are half-integers, that is {Qi} = {Q2} = ^ 


instead of (62). The leading order is still given by (70), specifically 


21 


not exponentially small, it can be read from the term j = 0 of (63) 


p =_i 
4 


1 


L-Lt -L, 


+ 0(e-2^) 


(71) 


with 


Lm — 1 ~ +- + !)+ 


H _ Qm 2 
2 Qm + h 


lQm\ 

lQm\ + 1 ’ 


m = 1, 2 . 


(72) 


We find here the generalization of the charge-symmetric case {Qi = Q2 = P P discussed in section^ Each charge 
qi and q 2 is screened by [QiJ and [Q 2 \ ions. The remaining counterion is free to roam in a region of size L — Li — 1 / 2 , 
and with zero electric field. This ion contributes to the pressure with a term . Here Li = (^lqij)oo is the 

size of the screening layer of [QiJ counterions formed around qi and L 2 = limL^oo(A — (x 7 v+i-[Q 2 J )) fhe size of the 
layer of [Q 2 \ counterions formed around q 2 (compare ( [T^ to (§, when [QiJ = [Q 2 \ =p). 

Like-charges that are natural integer multiples of — e. In this case, the screening is not frustrated as in 
all the previous situations. Simply Qi counterions will screen the charge qi forming a neutral object, and similarly 
around q 2 there will be a screening cloud of Q 2 counterions. Since both objects with their screening clouds are neutral, 
the effective force between them is expected to be Pc ^ 0^. This can be verified from the expression for the partition 
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function applicable here, Eq. ([64|). If L -> oo, we have 

Qi Q 2 


-L 


Zr = 


1 4e 


2 Q 2 + 1 

2 VQi(Qi + 1) Q2{Q2 + 1) 


1 ^ 2gi + i 

-L/ ~r ~ 


- 1 


+ 0(e-"^) 


Therefore, 


Pr = 4e 


-L Ql Q2 

Qi + 1 ^2 T 1 


^ ^ 1 ^ 2^1 + 1 _2^2 + 1 _ 2 


2 \gi(gi +1) Q2{Q2 +1) 


+ 0(e-^^). 


(73) 


(74) 


D. Density profile 


With the above results, we can obtain an explicit expression for the density profile of counterions 


r 

‘ •^k —1 


n{x) = 


<Xk-l<Xk=X<Xk + l<---<XN 




Zc{N^ L, Qi, Q 2 ) 


(75) 


Notice that due to the fact that each particle only feels a constant electric field proportional to the difference between 
the number of charges which are at its left and right sides, the potential energy has the following property 


U{N,L,Qi,Q 2) = U{k -l,Xk,Qi,Q2- {N -k + l)) + U{N -k,L - Xk,Qi- k,Q2 ). 


(76) 


This can be interpreted as follows. If the particle at position Xk is fixed, the system decouples into two independent 
systems, one of size Xk with k — 1 particles, and the other one of size L — x^ with N — k particles, with the appropriate 
charges at each boundary (obtained by summing the charges at the left side and right sides of Xk of the original 
system). Then, the computation of the integrals in (75) simply yields the product of the two partition functions of 
each subsystem. 


Ylk=i - 1 ,x,Qi,Q 2 - N ^ k - l)Zc{N - k,L - x,Qi - k, Q 2 ) 

Zc{N, L, Qi,Q2) 


n{x) = 


where each Zc should be replaced by its appropriate corresponding expression from (62) or (68). 


(77) 


1. Contact density and pressure 

From this expression we can verify the known relation between the contact density at x = 0 (or x = L) and the 
pressure m- Indeed, notice that 


^ Zc{N - 1, L, Qi - 1, Q 2 ) 

Z,iN,L,Qi,Q2) 


(78) 


On the other hand, from Eq. (62) we can verify that 


dZc{N^ L, Qi, Q 2 ) 


Zc{N — I, L, Qi — I, Q 2 ) — (Qi)^ , 


(79) 


where this last relation was obtained by writing —{j — = {N — j){j + a) 

that Qi = {N a)/2. Therefore, we find 


((A^ + a)/2)^ in (62), and recalling 


Pc = n(0) - (gi)2 = n(L) - {Q^f . 


( 80 ) 


The last equality is obtained using the same argument on x = L in n(x). 
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FIG. 8: The density profile for N = 2b and N = 26 counter-ions and L = 10. Notice that in the case where the number of 
counter-ions is odd, N = 25, the density far from the edges converges to a non zero value 1/(L — 2(x)oo), here close to 0.124. 


2. Asymptotic behavior of the density 


Let us consider the case a = 0, ie. Qi = Q 2 = N/2. Figure shows a plot of the density profile for N = 2b and 
N = 26. Notice that in the case N = 26 even, the density falls off quickly to zero far from the boundaries x = 0 
and X = L. On the other hand, when = 25 is odd, the density does not fall to zero, but goes to a non-vanishing 
value shown by the horizontal line. This corresponds to the density of the free counterion, responsible for the effective 
attraction between the two charges qi and q 2 as discussed earlier. 

To quantify this behavior, consider expression (77) for the density in the case = 2p + 1, and Qi = Q 2 = P + 


n{x) = 


-l,x,p+l,k-p- |)Zc(2p+ 1 - k,L - x,p - k ^ l,p^ 1) 

^c(2p + l,I/,p+|,p+^) 


(81) 


In this sum, the partition function Zc{k — 1, /c — p — |) corresponds to a system with charges —e{p + and 

—e{k—p— |) at its boundaries. If /c < p, these two charges carry opposite signs, therefore, Zc(/c —l,x,p+|,/c—p — |) 
is given by Eq. (62). Then, if 1 ^ x ^ L, Zc{k — l,x,p+ — p — |) = On the other hand, the 


second partition function, Zc{2p-\-l — k^L — x^p — k \ ^p A- ^), corresponds to a system with charges —e{p — k 
and —e{p + |) at its edges, li k < p^ these two charges carry the same sign and are half integers multiples of e, 
therefore Zc{2p + 1 — k^L — x^p — k A- A- \) should be obtained by using Eq. (64). In particular one can notice 
that if 1 <C X ^ L, then Zc{2p + 1 — /c, L — x^p — k A- \^p + |) = Therefore, in the sum ([sij) all terms 

with k <p decay exponentially fast when x is far from the boundaries: they are of order O 


same argument could be applied to all the terms with /c > p + 2, with the roles of Zc{k — 1, /c — p — |) and 


Zc{2p +1 — k^L — x^p — k A- \ ^pA- \) interchanged. Then, only one term in the sum (81) survives, it corresponds to 


k = p+1, which is precisely the index of the position of the free counterion. In this term, both Zc{k — 1, k—p—^) 

and Zc{2p +1 — — x,p — with /c = p + 1, correspond to a system with charges —e{p + and e/2 at 

its edges (notice the opposite signs), and those partition functions should both be computed using ( [6^ . The leading 
order of these partition functions, when 1 <C x ^ L, is 


/ 1 lx 

Zc(p,x,p+ -,--) 


^-xl4 2 2 ^-{L-x)l4 

,/ , -.x, and Zc(p,L-x,--,p+-) ^ 

p!(p+l)! 2 2 p!(p+l)! 


( 82 ) 
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while the leading order of the denominator of (81) is 


1 1 

Z,(2 p+l,L,p+-,p+-)~ (^!(^^1)!)2 


L-2- 


p +1 


(83) 


This gives 


n{x) 




for 1 <C X <C 1/. 


(84) 


This is the analytical confirmation of the intuitive analysis of section IIA where it was explained that when N is 
odd, there is one free ion roaming between the two charges with an available space equal to L — 2{xp)oo, as shown in 
figure 1^ 

In the case where N is even, a similar analysis shows that all terms of the sum ([^ fall of exponentially fast when 
X is far from the boundaries. 


E. The large N limit 

It is interesting to consider the limit N oo. Due to the electroneutrality condition ^ 1 +^ 2 + eA" = 0, one needs 
to consider different situations: whether qi and q 2 are kept finite, then necessarily the charge of the counterions e 
should vanish as 1/A. Then we notice that this is also a mean field regime. The other possible limit is to consider 
that e has a non vanishing finite value, then qi and/or q 2 should go to infinity as A. 


1. Mean field limit, A ^ 00 and e ^ 0. 


Momentarily, it is best to return to dimensional units L and P: the rescaling by is not appropriate here, because 
e ^ 0. Consider the equation of state (57) derived in the isobaric ensemble, which now reads 


N 


m = Yi 


1 


^(e^ + g<)2+P 


i=o 


- f 

e Ja. 


dy 


(85) 


where </< = —e(3< and = —eQ>. Since e —> 0, one can recognize a Riemann sum and replace it by an integral. 
This finally leads to 


92 


pe{L)Vp = arctan + arctan . 

VP VP 


( 86 ) 


We recover here the implicit relation between (L) and P from the mean field theory as described by the Poisson- 
Boltzmann equation Indeed, referring for instance to where the mean field regime of the present problem 

was considered, Eq. (86) can be directly obtained from a simple linear combination of Eqs. (16) and (17) of [27]. Notice 
that the interesting effects, such as like-charge attraction, stemming from the discrete nature of the charges, are lost 
in this mean field limit. Like-charges will always have a repulsive effective interaction in the mean field regime [2IH23]. 
A related comment is that the asymptotic negative pressure reported for odd A in section]^ P = —g^/A^, vanishes 
in the limiting process addressed here. 

It should be noted that the present limit is also the thermodynamic limit, since we have to remember that e is of 
order 1/A, therefore in the left hand side of (86) (P) should be of order A. To make this more apparent, introduce 
the average distance per ion {£) = {L)/N (inverse of the density), then (86) becomes 


P{Q1^Q2 


P = arctan arctan . 

VP VP 


(87) 


2. Limit A ^ oo and e fixed. 


In this situation, the charges at the edges qi and q 2 should be of order A, or at least one of them. Consider the 
case when both Qi > 0 and Q 2 > 0 are of order A. Then, when A ^ 00 , Eq. (57) can be put in the following form 
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by shifting the index of the summation by [Q<J, 


oo ^ 

(e - so ^ 


£= — oo 


{i-{Q<}r + P' 


( 88 ) 


Notice that by shifting the index £ by one, we can replace {Q<} by {Q>} if necessary. One can then write 

1 


= E TTTTS 


{£ - min({Qi}, {< 32 }))^ + P 


(89) 


Notice that in this analysis, the limit depends on how Qi and Q 2 are taken to infinity, and assumes that the fractional 
part of them is kept fixed as N is increased. 

To cover the whole range of values for (L) from 0 to +oc, it is necessary that P covers the range from 
— min({(5i}, {( 52 })^ to + 00 . We recover the same phenomenology as in the case N finite, when (L) ^ 00 , 
P — min({(5i}, {( 52 })^- So, the pressure can become attractive, except in the case where Qi and Q 2 are inte¬ 
gers. Eq. ([8^ can be made more explicit in two particular cases. When Qi and Q 2 are integers. 


{L) 



1 

'pTp 


TT coth(7r\/P) 


and when Qi and Q 2 are half integers. 


(L) 



1 

(TTIFTp 


TT tanh(7r\/P) 


(90) 


(91) 


When Qi and Q 2 are not integers, the value of (L) for which the pressure changes of sign is given by putting P = 0 




1 

(^-min({Qi},{Q2}))^ 


— + '^^({^ 2 }) • 


When Qi and Q 2 are half-integers this reduces to (L*) = tt^. 


(92) 


IV. CONCLUSION 

We have studied a simple one-dimensional system as a model to understand the effective interaction between 
charged particles that are screened by counterions only. This model evidences the possibility of attraction between 
two like-charges at large separation. The physical phenomenon behind this attraction is a frustration of the screening 
process due to the discrete nature of the electric charges. More specifically, if the two like-charges are not integers 
multiples of the charge of the counterions, a perfect screening of the charges is not possible, and there will be a 
“misfit” counterion, responsible for the over-screening of one of the like-charges, leading to an effective attractive 
force. A by-product is that in the mean-field limit where discreteness effects are washed out, no like-charge attraction 
is possible, a well-known phenomenon. 

The present model is in addition interesting from a purely theoretical perspective, since it is exactly solvable: it 
is possible to compute explicitly its partition functions (isobaric and canonical), the pressure (effective force) and 
the density profile of the counterions. Although the specific exact results and expression for the effective force are 
particular to this one-dimensional model, the physical mechanism responsible for the attraction between like-charges 
could also be applicable for three dimensional situations [19]. In particular the case N = 1 leads to an equation of 
state that is equivalent to that found under strong coupling for three dimensional planar interfaces, screened by point 
counter-ions interacting through the standard 1/r Coulomb potential [TTl[T^[T7] . 
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Appendix A: Two equal charges: canonical expressions 


The inverse Laplace transform can be computed with integral inversion formula which can be evaluated using the 
residue theorem 




yk = 0 + i) + P 

Each residue is straightforward to compute 


w=E R?i, 




TIP 


llk=0 (^+1) 


1 2 


^PL 


Res 




^ nfe=o (^+5) Y{k=o,k^j i^+l) “ 0 + 5) 

Writing 


z=^j (^ + ^) - {j + I) 


1 


1 


1 ^ 1 


1 


(A:+- (j + 1)2 ^k-j){k^j^l) 2j^l\k-j k^j^lj' 

the above product and sum can be simplified 

1 _ i-iyi2j + i) 


Ylk=o,kjij (^ + 5 ) “ 0 + 5 ) 


{p-j)Kp + j + 1)! ’ 


and 


E 


1 _ 2 f 1 i\ 

+i)' - (i+i)' “ 1 1 ^) 


(Al) 


(A2) 


(A3) 


(A4) 


(A5) 


Gathering all results, the exact explicit result for the canonical partition function is found in the form of Eq. (36). 
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